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This  technical  note  describes  the  development  of  a computer  program 
for  the  calculation  of  multi-ping  probability  of  detection  as  a part  of  a 
total  computer  model  for  sonar  prediction.  This  memorandum  has  been 
prepared  because  it  is  believed  that  the  information  may  be  useful  in  this 
form  to  others  at  Naval  Undersea  Research  and  Development  Center.  It 
should  not  be  construed  as  a report  as  its  only  function  is  to  present  a 
portion  of  the  work  as  information  for  others.  This  technical  note 
supersedes  NUWC  TN  90. 

This  work  has  been  supported  by  NAVSHIPS  Code  00V1  under  subproject 


SF  11-121-105  Task  14318. 


INTRODUCTION 


v/ 


One  type  of  detection  criterion  of  interest  for  use  in  sonar  systems 
is  a "counting  criterion"  where  multiple  looks  at  a target  area  are  avail- 
able. The  receiver  output  for  each  look  or  ping  is  compared  with  a thresh- 
old and  generates  a binary  one  output  if  that  threshold  is  exceeded  or  a 
binary  zero  if  it  fails  to  exceed  the  threshold.  The  multiple  looks  will 
now  appear  as  a binary  sequence  and  the  "counting  criterion"  examines  a 
group  of  length  m,  counting  the  number  of  ones  in  that  group.  If  the 
number  of  ones,  or  threshold  crossings,  equals  or  exceeds  the  integer  n 
the  detection  criterion  is  satisfied  and  a target  is  declared  present. 

Otherwise  another  ping  in  the  sequence  is  examined.  This  counting 
criterion  is  often  called  an  n out  of  m criterion.  Since  successive 
groups  of  the  m pings  have  m-1  pings  in  common  the  probability  of 
satisfying  the  n of  m criterion  at  any  'given  point  of  the  sequence  is 
dependent  upon  the  m-1  past  outcomes.  Such  a sequence  is  often  called 
a Markov  chain  of  m— 1 order.  It  is  of  interest  to  examine  the  joint 
probability  of  such  overlapping  groups  of  independent  pings  and  to  determine 
the  probability  that  n or  more  of  the  group  (in  any  order)  exceed  a fixed 
threshold,  thus  qualifying  that  group  as  satisfying  an  n out  of  m detection 
criterion.  In  addition,  the  cumulative  probability  that  an£  group  in  the 
sequence  satisfies  the  n of  m criterion  is  desired.  Since  adjacent 
groups  have  many  pings  in  common  the  groups  themselves  are  not  independent, 
even  though  the  pings  are  assumed  to  be.  It  is  assumed  that  the  single  ping 
probability  of  exceeding  the  threshold,  P^,  is  available  and  all  solutions 
will  be  in  terms  of  the  single  ping  detection  probability. 
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Figure  1:  Binary  Output  of  Threshold  Receiver 

t* 

JOINT  PROBABILITY  FOR  K GROUP 

Z(t)  is  a binary  function  which  represents  the  sequence  of  single 

ping  outputs.  If  the  kth  ping  which  is  received  at  time  tfc  exceeds 

a fixed  receiver  threshold,  Z(t^)  is  equal  to  one;  otherwise,  Z(t^)  is 

zero.  P,  is  the  single  ping  probability  that  the  k pulse  exceeds  the 

threshold  and  thus  the  probability  that  Z(tk>  will  equal  one.  We  now 

s t 

examine  a group  of  m adjacent  pulses  consisting  of  the  k-m+1  to 
the  kth  and  write  down  the  probability  that  Z(tk)  for  this  group  will 
contain  exactly  n ones  and  m— n zeros  under  the  assumption  of 
independence  between  pings. 

If  the  P.'s  were  equal,  the  expression  for  the  probability  of  exactly 
K. 

n out  of  m is  the  familiar  binomial  form: 

P{n/m,  k)  = Pn(l-P)m^n,  P = Pk  for  all  k 

I « 

and  the  probability  of  n or  more  out  of  m is 


P{>  n/m,  k)  - f ^Pj(l-P)mj 
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For  the  more  general  case  of  unequal  P^'s  fP{n/m,k}  is  the  sum  of  all 
possible  products  of  n different  Pk's  and  m-n  different  (1-Pk)'s. 
Thus,  if  n = 3 , m = 4 . 

P{3/4,k}  = Pk  Pk_1  PR_2  d-pk_3)  + Pk  Pk-1  (1_Pk-2)  Pk-3 
+ Pk  (1-P^)  Pk_2  Pk_3  + (l~Pk)  Pk_x  pk_2  Pk-3* 

The  number  of  terms  in  the  sum  is 


(“)  A — 
\n  / * n 


m! 


! (m-n) ! 


In  the  above  example  P{^  3/4, k)  is  given  by 


P{>  3/4, k}  - P{3/4,k}  + P{4/4,k} 


Pk  pk_]^  pk_2  U pk_3^  pk  pk_^  ^-Pk-2^  Pk-3 


+ Pk  (1-P^)  Pk_2  Pk_3  + (1-Pk)  pk_x  pk_2  pk_3 


+ P P P P 
rk  k-1  k-2  k-3 


If  P = P the  above  expression  reduces  to 

K 
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P{^  3/4, k}  = P3  (1-P)  + P3  (1-P)  + P3  (1-P)  +P3  (1-P)  + P4 

= 4 P3  (1-P)  + P4 

“ (a)  p3  (1-P)  +(4)^ 

the  same  as  previously  noted. 

Cumulative  Probability: 

It  is  of  interest  to  calculate  the  probability  that  a string  of  pings 

of  length  k contains  a sub-string  of  length  m which  satisfies  the  n of 

m criteria.  This  sub-string  can  be  located  anywhere  in  the  string.  This 

probability  will  be  denoted  by  P ^ an<^  a nionotonically  increasing 

function  of  the  string  length  k.  To  obtain  a general  algorithm  for  the 

calculation  of  P , one  approach  is  to  generate  all  possible  binary 
cum 

functions,  Z(t)  of  length  k,  selecting  those  that  meet  the  criterion  and 

calculate  their  probability  and  sum.  Since  the  number  of  possible  sequences 

is  2 this  procedure  is  not  practical  for  large  k.  A second  approach 

in  which  P (k)  is  a function  of  P (k-1)  and  a smaller  set  of  binary 
cum  cum 

functions  will  be  described. 

First  a few  definitions.  Let  E(k)  be  the  event  that  a given  binary 

1 

function  of  length  k has  n or  more  ones  (successes)  occurring  within 

i ___ 

the  last  m positions  and  E(k)  the  complement  of  that  event.  Thus,  if. 

Z(t)  is  written  as  a binary  number  with  the  right  most  digit  representing 

Z(t^)  then  for  k = 7,  n=  2,  m=  3,  the  string  0010010  fails  while  0010101 

is  in  event  E(7).  With  this  notation  P (k)  is  the  probability  that 

cum 


- 1 (i) 


Fk  (1-P) 


k=3 
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E(k)  occurs  or  E(k-l)  or  E(k-2)  ...  or  E(m).  Using  DeMorgan's  law 


P (k)  = P(E(k)UE(k-l)UE(k-2)U- • -UE(m)) 
cum 

= P (E(k-l)UE(k-2)U* • •UE(m)) 


4-  P(E(k)nE(k-l)nE(k-2)n- • -E(m)) 

P (k)  = P (k-1)  + P(E(k),  E(k-l) , E(k-2) , • •■  , &(m)). 
cum  cum 

Using  the  notation  P(A/B)  for  the  probability  of  A given  that  B occurred, 

P (k)  is  given  by 
cum 

P (k)  = P (k-1)  + P(E(k) I E(k-l)  E(k-2) , . . . E(m) ) 
i cum  cum 

- P(E(k-l) , ¥(k-J)  , . . . E(m) ) 

but 


P(E(k-l) , E(k-2) , . . E(m) ) = P(E(k-l)UE(k-2)U. .UE(m) ) 

= 1-P  (k-1)  giving 

cum 

P (k)  = P (k-1)  + P (E(k)  | EOk-T)  , E(k-2), E(m))  ■ (1-P  (k-1)). 

cum  cum  1 cum 


*U  is  union,  n is  intersection  hereafter  written  with  a comma 
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The  E(k)  and  E(k-m),  have  no  common  pings  which  determine  the  events 
success  or  failure;  therefore. 


P(E(k) | E(k-l) , E(k-2) E(m))  = P (E(k) | E(k-l) , E (k-m+1)). 


Since  the  event  E(k-m+l)  depends  on  pings  number  (k-m+1)  back  to  (k-m+1) 

- m+1  or  (k-2m+2)  and  the  event  E(k)  depends  only  on  pings  number  k back 
to  k-m+1  the  P(E(k) |E(k-l),  E(k-l-m+l))  may  be  determined  by  examining 

all  possible  binary  sequences  of  length  k - (k-2m+2)  + 1 or  2m-l  ending 
with  the  k*"*1.  There  are  2^m  ^ possible  sequences  which  can  be  partitioned 
into  the  following  groups  using  n=2,  m=3  as  an  example. 


6 


m 


Only  those  strings  (of  length  2m*“2)  in  Group  I can  result  in  the  event  E(k)  |E(k  1) 

...  E(k-m-fl)  no  matter  what  the  outcome  of  the  k ping.  If  the  k ping  is 

a success  (1)  the  event  is  satisfied  by  the  strings  in  Group  I and  none  of  the 

th 

strings  in  Group  II  or  Group  III.  If  the  outcome  of  the  k ping  is  a failure 

r 
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(0)  the  event  fails  for  all  groups.  Therefore,  all  successful  strings  of 
length  2m- 1 are  obtained  by  taking  those  ip  Group  I and  adding  a one  in  the 

fch 

k position  in  the  example: 


The  probability  of  the  event  E(k)|E(k-l),  . . E(k“irH-l)  may  now  be  written 
in  terms  of  the  probability  of  the  strings  in  Group  I and  Group  II. 


P(E(k)|E(k-l), 


-p(E(k)-  E<k'1)’  E(k-gl+m 

P (E(k-l) , E(k-nri-l) ) 


p(k)  • (Sum  of  the  probability  of  all  strings  in  Group  I)* * 

(Sum  of  the  probability  of  all  strings  in  Group  I or  Group  II) 


P(k)  * P(Gp  I) 

P (Gp  I)  + P (Gp  II) 

th 

*Note:  In  all  successful  strings  the  k ping  must  be  a 1 (success)  and 

the  k~m  ping  a 0 (failure). 
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In  the  example  of  n=2,  m=3,  k=5,  and  P(i)  the  probability  that  the  i 
ping  exceeds  a fixed  threshold. 


th 


Pk(GP  I)  = (l-P(l))  • (1-P (2) ) • P (3)  • (1-P (4) ) 
+ P(l)  • (1-P (2) ) • (1— P (3 ) ) • P (4) 

+ (l-P(l))  * (1— P (2 ) ) • (1-P(3))  * P(4) 


Pk(GP  II)  = P(l)  • (1-P (2) ) • (1-P (3))  • (1-P(4) ) 

+ (l-P(l))  ' P (2)  • (1-P(3))  • (1-P (4) ) 

+ (l-P(l))  • (1-P (2) ) • (1-P(3))  • (1-P(4)) 


and 


P (k)  = P 
cum  cum 


(k-1)  + 


P(k)  * P,  (GP  I) 

{1-P  (k-1)} 

P (GP  I)  + P (GP  II)  CUm 


After  the  strings  of  length  2m- 2 in  Group  I and  Group  II  have  been 
generated  (for  a fixed  n and  m)  they  are  retained  for  calculating 
P(E(k))  |E(k-l)  , E(k-nrfl))  for  all  k > 2 m-1  by  changing  only  the 

individual  ping  probabilities,  P(k),  which  enter  into  the  calculation. 

A general  flow  diagram  for  the  calculation  of  Pcum(k)  f°r  a given  n 
and  m is  given  in  the  appendix.  The  writing  and  check-out  of  this 
program  was  completed  by  Mr.  K.  A.  Faucher.of  Code  556. 
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Appendix:  Flow  Diagram 

MAIN  PROGRAM 


cum 


START 


Input:  I 

M,  L 

P(k)  k-1, 2,  ..L 

K=0,  Q=( 

) 

is  an  array  of  length  L 

containing  the  cumulative 
probability, 
is  an  array  of  length  L 
containing  the  ping  • 
probabilities . 
is  an  array  of  length  Q 

containing  the  sequences 
in  Group  I. 

is  an  array  of  length  V 

containing  the  sequences 
in  Group  II. 


Call  PGEN  to 
calculate  the 
prob  {Gp  1} 


Call  PGEN  to 
calculate  the 
prob  {Gp  II} 


P (k)  =P  (k-1) 
cum  cum 

+(1-P  (k-1)) 

cum 

*P(GP  I)/P 

(GP  I)  + P (GP  II) 
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SUBROUTINE  INITIALIZE 


Initialize  generates  the 
sequences  in  Group  I and 


A 
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Q = 2 * Q 
q = 0 
sum  = 0 


3 


> 
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SUBROUTINE  PGEN 


PGEN  accepts  an  array  sequence 
(either  S or  T)  and  calculates 
the  sum  of  the  probability  of 
all  the  sequences  in  the  array 
using  P(k) . 


> 
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